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Abstract 

The quintessential vortex solution in (3+l)-dimensional nonrelativistic, incompressible fluid 
mechanics is the Burgers vortex. We show that, in a finite domain, conformal fluids also admit 
hot vortex solutions with everywhere nonrelativistic speeds. These are identical to Burgers' 
solution, except that their radius is reduced by a factor of 2/-\/3. A rough calculation indicates 
that at RHIC these vortices are indeed smaller than the fireball itself during thermalization. 
Similarly to the Burgers vortex, these solutions manifest vortex stretching which avoids short 
distance singularities and so suggests that conformal fluid flows share the same nonsingularity 
as solutions of the Navier-Stokes equations. Naively generalizing this calculation to an arbitrary 
equation of state w, we observe that the Burgers vortex radius diverges as w crosses —1. While it 
has been argued that such a crossover leads to an instability in certain perfect fluids, the absence 
of Burgers vortices and therefore vortex stretching suggests that, in addition to the well-studied 
big rip singularities, viscous phantom fluids generically develop vorticity singularities. 
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1 Motivation 



The initial conditions under which the evolution of the Navier-Stokes equations of fluid dy- 
namics lead to a nonsingular solution are unknown. However, as has been argued in Refs. [H 12], 
numerical simulations indicate that such singularities are avoided by a mechanism called vortex 
stretching. In particular, in contrast with Kolmogorov's scale-invariant model of turbulence [3], 
turbulent configurations contain random distributions of thin vortices These vortices are 
all deformations of an exact solution of (3+l)-dimensional incompressible hydrodynamics, the 
Burgers vortex [5]. As we will review below, the Burgers vortex constantly stretches, pumping 
vorticity out of potentially singular regions and thus cutting off velocity gradients, ultimately 
preserving the continuity of the flow. 

In this note we will extend Burgers' incompressible fluid solution to the case of a (3+1)- 
dimensional conformal fluid. We will demonstrate that, to leading order in an expansion in 
the inverse temperature, hot vortex solutions with nonrelativistic velocities exist, but with a 
radius which is smaller than that of the Burgers vortex by a factor of 2/ a/3. In the case of 
incompressible fluids, the fact that the constant stretching of Burgers vortices yields a nonsin- 
gular evolution of the flow has not been proven, but has been supported by numerous numerical 
simulations. In the conformal case, even such weak support for nonsingular evolution is absent. 
The statistical mechanics of webs of thin vortices is highly dependent upon the theory of their 
nonaxially-symmetric perturbations [1] which will not be discussed in this note, and so may in 
principle be quite different in the conformal case. However the fact that such similar solutions 
also exist and indeed stretch in the conformal case nonetheless implies that they do indeed carry 
vorticity out of high-vorticity regions, and so may also lead to the nonsingularity of conformal 
fluid flows. 

Such conformal fluids are of interest for a number of reasons. One reason is that, at temper- 
atures far above the deconfinement transition, the quark gluon plasma is a conformal fluid. It 
is a charged conformal fluid with various electromagnetic and color fields, however in heavy ion 
collision data analysis they are often [^j approximated by uncharged, conformal fluid£]. Near 
the deconfinement transition the conformal anomaly leads to a large bulk viscosity which de- 
stroys the conformal approximation, and so this approximation may be trusted at best during 
the initial evolution of man-made and cosmological quark gluon plasma. The rapid thermal- 

"'^More precisely, the baryon number density does not appear in the equations for the velocity and the 
thermodynamic variables, and so baryon number conservation may be imposed separately, while the effects of 
electric charge conservation are negligible [7]. 
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ization of the quark gluon plasma at heavy ion colhsions suggests a period of turbulence, and 
the references above suggest that, at least in nonrelativistic (3+l)-dimensional turbulent flows. 
Burgers vortices play a crucial role in the redistribution of vorticity This thermalization occurs 
when the quark gluon plasma is much hotter than the deconfinement transition temperature, 
and so the conformal approximation can be trusted. 

Not only does quark gluon plasma appear on Long Island and near Geneva, but also our 
universe was once made of quark gluon plasma. Both heavy ion and cosmological data so far 
are consistent with the perfect fluid approximation, whereas viscosity plays a crucial role in the 
Burgers vortex. However the nature of the turbulent flow plays a key role in understanding the 
evolution of the initial data of heavy ion collisions to the quark-gluon plasma phase, which is 
perhaps the greatest roadblock to data analysis. Thus one may hope that data on the structure 
of the vorticities is not entirely out of reach. 

Another reason for interest in such conformal fluids is that it has recently been demonstrated 
[8] that there is a one-to-one correspondence between their flows and hot black brane solutions 
to Einstein gravity with a negative cosmological constant. The vortex solution that will be 
found below lies within the regime of validity of this correspondence, and so is guaranteed to 
provide a new solution to Einstein gravity [H] . However, as we will consider the fluid flow only 
in a finite region of spacetime, the gravity dual will not be geodesically complete. 

We begin in Sec. |5] with a review of the Burgers vortex solution of the Navier-Stokes equa- 
tions for an incompressible fluid. In Sec. [3] we generalize this result to hot conformal fluids, 
solving the relativistic Navier-Stokes equations to leading order in a v/c expansion. We find 
that indeed there exist, in a local domain, solutions analogous to the Burgers vortex with a 
velocity that is everywhere nonrelativistic. These solutions are characterized by their vorticity 
profiles, which we see are identical to those of the Burgers vortex but spatially reduced by a 
factor of We also determine the spatial profiles of the various thermodynamic quanti- 

ties, but stress that these may receive large corrections from relativistic effects. Finally Sec. H] 
contains more speculative topics. Ignoring the background expansion, we use crude estimates 
of parameters of the RHIC fireball to determine the size of a potential Burgers vortex. This 
rough estimate suggests that such a vortex indeed would fit within the fireball, however it 
represents only a slight improvement upon a straight dimensional analysis. We also generalize 
this construction to an arbitrary equation of state, finding that phantom fluids do not possess 
Burgers vortices, and so cannot benefit from the desingularization of vorticity caused by vortex 
stretching. 
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2 Review of Burgers' Solution 



We will start by reviewing the Burgers vortex solutions to the Navier-Stokes equations for 
a nonrelativistic, incompressible fluid. We will be interested in time-independent solutions in 
3 spatial dimensions, x, y and z, with a divergenceless background constant strain field. The 
solutions will have an axial symmetry which rotates the x — y plane, and the strain field will 
point inward along that plane with a magnitude g. In other words, there will be background 
flow which compresses the fluid along the x — y plane, and extends it in the z direction, with 
a velocity that scales linearly with distance. The flow is homogeneous in the sense that it 
is translation invariant up to a Galilean transformation. Clearly such a condition may not 
be imposed upon a relativistic system, as the corresponding Lorentz transformations would 
commute only up to a Thomas precession term, and so a strain field at relativistic speeds may 
at most be homogeneous either in the x — y plane or along the z direction. However, in the 
present note, we will restrict our attention to flows in which the velocity is always nonrelativistic. 
A generalization to the relativistic case should be straightforward, mirroring the analysis of the 
(2+l)-dimensional case in Ref. |lUj . 

The fact that the velocity varies linearly with respect to the distance in a constant strain field 
implies that at sufficient large distances, the velocity nonetheless becomes relativistic. Therefore 
we will restrict our attention to a bounded spatial region. Physically this is a reasonable 
restriction when studying turbulent flows, because the strain fields are only approximately 
constant in local domains. While the individual particle velocities in cosmological and heavy 
ion quark gluon plasmas are heavily relativistic, it is not at present clear whether the flow 
velocities are also relativistic. In addition, both cosmological and man-made quark gluon 
plasmas expand rapidly. Our divergence-free velocity approximation is only valid if the shear 
g locally dominates this overall expansion. For example, in the case of the RHIC plasma, this 
condition is g » 10^^ /sec fTT]. In fact, without such a large shear, the vortex radius that we 
will find would barely fit inside of the quark gluon plasma fireball. 

The axially-symmetric Burgers vortex solution consists of a superposition of the strain 
field, described in the previous paragraph, with a rotational flow with angular velocity /, 
which is a function of the cylindrical radial coordinate p = a/x^ + y'^. This function breaks the 
homogeneity of the flow, leaving a U{1) rotational symmetry as well as translation symmetries 
in the z and time directions. The total velocity v(/9) is the sum of the strain field and the 
rotational flow 

Vx{p) = -xg - yf{p), Vy{p) = xf{p) - yg, v^{p) = 2gz. (2.1) 
The time- independent Navier-Stokes equation consists of 3 components indexed by the coordi- 
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nate label k 

= vWjv'' + d^V{p, z) - ud^djv'' (2.2) 

where V{p,z) is the pressure per unit density and u is the kinematic viscosity, which is taken 
to be constant as in this section we are interested in Newtonian fluids. 

The cylindrical symmetry implies that it is suflicient to consider only the x and z components 
of the Navier-Stokes equations (12. 2p . The z component is quite simple, as the linearity Vz implies 
that its' second derivatives vanish. It is simply 

dzV{p,z) = -v'd,v, = -4gh (2.3) 

which is solved by 

V{p,z)=Vo{p)-2g^z'. (2.4) 

In other words, the pressure does not respect the z translational symmetry of the system. Total 
shifts of the pressure are irrelevant in incompressible fluids, however this implies that in the 
relativistic case the z-translational symmetry will necessarily be broken unless the p-dependence 
of Vz and Vp can be chosen so that this decrease in pressure is precisely compensated by a 
relativistic 7 factor. If there is no such compensation then, in the case of relativistic conformal 
fluids, the kinematic viscosity will be minimized in the middle of the vortex, and so the Burgers 
vortex radius will have a minimum. 

The only remaining component of the Navier-Stokes equation (12. 2p is the x component, 
which is 

= x[-fipr + 9' + V',ip)/p] + y[2gfip) + pig + 3u/p')fip) + ufip)]. (2.5) 

For each value of p, this must hold everywhere on the corresponding circle of the x—y plane, and 
so both quantities in the square brackets must vanish separately. This leads to two equations. 
The first equation may be integrated to yield the pressure 'Po(p) as a function of a constant of 
integration Vqo, the angular velocity /(p) and the strain g 

Vo{p)=Voo+ f dp'p'{f{p')-g^). (2.6) 
Jo 

In particular, at large p we will see that the angular velocity tends to zero, and so the pressure 
per density tends to Pqo ~ The nonrelativistic approximation is violated before the 

pressure per density difference is of order one, however the pressure per density is very small 
in a nonrelativistic, incompressible liquid. Thus one may expect this approximation to break 
down at radii such that the velocity is of order C( V^oo); well below the speed of light. 

The vanishing of the second square bracket in (12. 5p may be reexpressed in terms of the 
vorticity 

^(P) = -5p(pV(p)). (2.7) 
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This is easily integrated to yield the Burgers solution [5] 

u;{p) = ce"^ (2.8) 

where c is an arbitrary constant of integration. Note that at large radii p the Burgers solution 
asymptotically tends to a constant strain field. The vorticity is concentrated within a vortex 
of radius 

PB = VWg. (2.9) 

3 Vortices in Conformal Fluids 

In the case of a conformal fiuid, even at nonrelativistic velocities, the pressure's contribution 
to the moment of inertia is no longer negligible. Thus the system is not well described by the 
Navier-Stokes equation. However translation invariance implies that the stress-energy tensor 
is still conserved. We will consider fiuids with no other conserved charges, and so the Landau 
frame fiuid fiow will be entirely characterized by the conservation of the stress tensor [12] 

d^T"" = 0, T^"-" = P{7]^"' + Au^'u") - ^a^r (3.1) 

Here P is the pressur^, m = f 7 is the relativistic 4- velocity, r/ is the dynamic shear viscosity, and 
a^^ is the shear. The first term in the stress tensor describes the perfect fiuid approximation, 
while the second incorporates viscosity. All functional dependencies have been suppressed, but 
they are analogous to those of the incompressible case, and we have used the conformality of 
the system to replace the density with three times the pressure. To generalize the Ansatz (12.1 p 
to this setting, we will define the constant g and function / by choosing the 4-velocity 

Ux = -xg - yf, Uy = xf- yg, = 2gz, m° = -uq = 7 = ^/l + p'^{p + g'^) +Ag^z^. 

(3.2) 

Note the divergence of the 4-velocity is zero and so the shear may be expressed as 

a^"" = P^'P{d,u^ + d^Up)P'"', P>"' = rj^"" + u'^u" . (3.3) 

In this note we will be interested in conformal fiows at nonrelativistic speeds. Therefore we 
will expand our results in v/c and only keep the leading orders. There are now 4 equations 
to consider, the 4 components of Eq. (13. ip . The three spatial components are the relativistic 
generalizations of the corresponding spatial components of Navier-Stokes, while the time com- 
ponent expresses the conservation of energy. The conservation of energy equation consists of 
the derivatives of the two terms in Eq. (13. ip . Assuming that the temperature and therefore 

is not the pressure per density V, which in the case of a conformal fluid is equal to 1/3. 
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the density and pressure are time- independent, at leading order in the velocity, the derivative 
of the first term is proportional to the divergence of u and so vanishes. The derivative of the 
second term does not vanish, however it contains elements of the strain with one temporal and 
one spatial index. Such mixings are suppressed by a factor of the velocity, and therefore this 
term is of order 0{v'^/c'^) and we will neglect it. 

In a conformal fluid every transport coefficient is equal to a dimensionless constant multi- 
plied by a power of the temperature T. We will make the standard assumption that all particle 
species are in local equilibrium, so that a single temperature is locally well defined. This ap- 
pears to be justified in the case quark-gluon plasma [13]. In particular, conformal fiuids are 
characterized by two dimensionless constants and cp such that 

P = CpT\ r] = c,T\ (3.4) 

In the case of the quark gluon plasma, RHIC data suggests Cr^ < Cp [H]. We will be interested 
in a very hot plasma, in which the temperature T is much greater than the scales of the various 
derivatives dv/v and dT/T. In particular this implies that the temperature is much greater 
than the shear T » g. In this case, the first term in the stress tensor (13. ip is much larger than 
the second. Indeed, without such an assumption, the derivative expansion of the velocity used 
in the form of the stress tensor (13. ip would be inapplicable, and the system would not behave 
like a fiuid. 

To leading order in the high temperature, low velocity expansion, the z component of the 
equation of motion (13. ip may be approximated by the derivative of only the perfect fiuid term 
of the stress tensor. Imposing time-independence, the conservation of 2;-momentum becomes 

= dk[P{r]''^ + 4uV)] ~ d,P + APdk{2gzu'') = d,P + IGg^zP (3.5) 

which is easily integrated to yield 

P = Pqc-^^'z^ (36) 

where Pq only depends on p. This result is superficially similar to the incompressible result 
(12. 4p . However one much recall that in this case the pressure per density V is constant, and 
so the density exhibits the same exponential falloff. Using the relations (13. 4p one may find the 
2-dependence of the temperature T and dynamic shear viscosity rj 

T = Toe-'^"^\ V = Voe-''"'" (3.7) 

where Tq and rjo only depend on p, and in fact themselves obey (13. 4|) along with Pq. These 
relations are of limited use because only the constant and conceivabljo the quadratic terms in 
(13.61) are reliable, the others receive relativistic corrections. 



■^More precisely, using the Ansatz p. 21) the quadratic terms are correct. However, as will be described in 
Sec. m a fully relativistic treatment requires an Ansatz with at least one additional degree of freedom, which 
has the potential to modify and even eliminate these quadratic corrections. 
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Again axial symmetry implies that the x and y equations of motion are equivalent, and so 
it will suffice to impose the conservation of x-momentum. In this equation the strain tensor 
(13.31) will enter at leading order. Up to cubic order in / and g, the spatial components of the 
strain tensor a defined in Eq. (13.31) are simply those of the nonrelativistic strain tensor du. For 
example 

a-- = -2g-^f' + 0(v^), a"^= ~ f + 0(v^), a'^' = 0(v^). (3.8) 
P P 

Therefore, in the nonrelativistic approximation, the conservation of x-momentum is the condi- 
tion 

= d^P + Au'dkiPu^) - d^iva^^) - dy{r^a^y). (3.9) 

The P'u^ term is subdominant to the P' term in the nonrelativistic regime, and so may be 
dropped. As was the case with the x component of the Navier-Stokes equation, this equation 
transforms as a vector under the SO (2) rotation symmetry, and so may be decomposed into 
a quantity multiplied by x and a quantity multiplied by y. Again, these two quantities must 
vanish independently. 

Using (13. 4p to reexpress the viscosity and pressure in terms of the temperature, the vanishing 
of the X coefficient yields 

= {AC,T^ + Qcr,T^g)T,' + AC,T^p{g^ - f) (3.10) 

where Tq' is the derivative of the temperature with respect to the cylindrical radial coordinate 
p. In the high temperature regime Tq » g this is easily integrated to yield a global solution 
for the temperature 

T = Tooe-2^'^'+/'^^''''(^'-^'). (3.11) 

As was the case with Eq. (13. 6p . this resembles the nonrelativistic result (12.61) for the pressure per 
density ratio. However, the pressure per density ratio V is now constant. Eq. (13. lip implies that 
at large distances in both the z and p directions, the temperature and therefore the viscosity 
and pressure shrink. Eventually this invalidates the high temperature approximation, but only 
beyond the nonrelativistic regime where (13. lip may be trusted. 

The final constraint comes from the y coefficient of (13. 9p . which is the local conservation 
of angular momentum or more precisely the condition that all of the sources of torque must 
cancel in a stationary configuration. Using (13.41) . in the nonrelativistic approximation it is 

= Ac^T\2fg + pgf) + 3c,/'TVp + c^T\f"T + 3/'T') (3.12) 

where all derivatives are taken with respect to p. This disagrees with its nonrelativistic coun- 
terpart (12.51) in two respects. First, the last term, which measures the spatial derivative of the 
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viscosity, is new. Using the p-dependence of T in (13. lip , one sees that this derivative yields a 
factor of p{p — g"^). This means that the f"T term is larger than the 3f'T' term by a factor 
of order the inverse velocity squared, and so the latter may be dropped in this nonrelativistic 
approximation. 

The second difference is that the kinematic viscosity u has been replaced with 

T% = ?^ (3-13) 

where we have used the fact that the density of a conformal fluid is three times its pressure to 
define a kinematic viscosity u. However the new factor of 3/4 in front of the viscosity changes 
the Burgers solution (12. 8p to 

uj{p) = ce"^ (3.14) 



and so the Burgers radius has decreased from the incompressible value of yvjg. This is the 
main result of our analysis, the radius of the Burgers vortex in a conformal fluid subjected to 
a constant strain g is 

1/3^ , , 

4 RHIC, Relativistic Vortices and Phantoms 

Using an initial RHIC dynamic viscosity r/ of order 10//m^, density of order 1 GeV/ fm? 
and strain of order 10^'^ /sec leads to a Burgers radius of ~ -1 — Ifm. The Burgers vortex 
radius is inversely proportional to the temperature, and so increases by an order of magnitude 
prior to freeze-out, eventually being about the same size as the entire fireball and invalidating 
the solution. Needless to say, the fact that the RHIC fireball expands means that the strain is 
not divergenceless, as has been assumed here. Burgers vortices in an expanding fluid will be 
treated elsewhere. Such a treatment is a prerequisite for a determination of the observational 
consequences (or lack thereof) of these vortices in the quark gluon plasma at the LHC. 

The RHIC quark gluon plasma achieves thermal equilibrium, and so one may define a 
(convention-dependent) local velocity. It is not known whether this velocity is relativistic, but 
if it is, then one must confront the inconsistency of our Ansatz (13.21) . In principle a consistent 
Ansatz requires that the coefficient functions in the velocity components possess a nontrivial 
z-dependence. However if the fully relativistic solutions share a ^-translation symmetry up 
to a Lorentz boost, generalizing the longitudinal symmetry of Burgers solutions, then the z- 
dependence will be dictated entirely by this symmetry and a constant characterizing the strain. 
Our Ansatz (13.21) is nonetheless inconsistent, as there are four equations of motion but only two 
free functions / and T of p. We have argued that to leading order in the velocity our Ansatz 
contains solutions, which are justed rescaled Burgers vortices. 
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In a fully relativistic treatment the velocity components Uz and Up must also be functions 
of p determined by the equations of motion. As is already a function of z, this naively leads 
to partial differential equations, which would be quite difficult to solve numerically. However, 
if one fixes at ^ = 0, then the value of Uz may be determined everywhere by imposing 
the aforementioned longitudinal translation/boost symmetry. This leaves four ordinary dif- 
ferential equations at 2; = and four undetermined functions of p: the longitudinal, angular 
and radial velocities and the temperature at z = 0. It therefore seems plausible that such a 
translation/boost-symmetric vortex exists. Note however that in this counting we have allowed 
M2 at 2; = to be an arbitrary function of p. Therefore it seems likely that this function will nec- 
essarily be nonzero, and so the relativistic vortex will not share the z ^ —z refiection symmetry 
of the incompressible Burgers vortex, instead these vortices would come with orientations. 

What would be the physical interpretation of Uz{p) at 2; = 0? Naively it seems like a center 
of mass fiow, however by shifting the coordinate z, one can arrange for the net z-fiow at z = 
to vanish, for example leaving T*^(z = 0) positive for p < Pc and negative for p > pc where pc 
depends entirely on the vorticity, background shear and background temperature. Of course, 
even if such bizarre solutions exist, it does not mean that they play an important role in the 
hydrodynamics of these fiuids. The above degree of freedom count merely suggests that such 
a feature may be necessary if one demands the same longitudinal translation/boost symmetry 
enjoyed by vortices in incompressible fiuids. 

One may readily generalize these vortices to a nonconformal fiuid with density e, no bulk 
viscosity and an equation of state 

P , . 

w = -. 4.1 

e 

We will not assume that w is a universal constant for a given fiuid, as negative values of w would 
then lead to an imaginary speed of sound and so an instability. However, motivated by the 
small fractional change in temperature seen above over the radial profile of the vortex solution, 
we will assume that w assumes a constant value on the entire vortex solution. Then, again 
dropping the sub dominant derivative of the viscosity, the conservation of angular momentum 
fl3.12p becomes 

= (^) P9{2f + pg') + v{— + /") (4.2) 
\ w J p 

and so the Burgers radius is 

In the incompressible and conformal cases w = and w = 1/3 this yields Eqs. (12. 9 p and ( I3.15P 
as it must. 

Not surprisingly, there is no vortex solution for dark energy w = —1. However Burgers 
vortices also do not exist for phantom fiuids {w < —1), unless the kinetic viscosity v is neg- 
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ative. In fact, without such a negative the conservation of angular momentum equation 
fl3.12p imphes that the viscosity increases, rather than decreasing, the angular velocity in such 
configurations and therefore probably leads to a singularity, or at least a breakdown of local 
equilibrium. Ignoring such potential instabilities, if u remains positive there will be no Burgers 
vortex solutions and so the vortex stretching mechanism for avoiding singularities in incom- 
pressible fluids does not apply to phantom fluids. Thus, it seems quite plausible that such 
fluid flows develop singular vorticities where the local equilibrium conditions and therefore the 
hydrodynamic description itself break down. 

Phantom fluids do not satisfy the null energy condition. A large class of microscopic theories 
described by such fluids have been shown to either be unstable or to admit superluminal 
propagatioiu in Ref. In the case of perfect fluids microscopically described by a single 

scalar fleld with a reasonably general Lagrangian, obstructions to crossinjf] w = —1 have 
already been described in Refs. [191 ED [22] • However the present argument is generic to any 
fluid with a positive kinematic shear viscosity and w < —1, even if it has never crossed w = —1 
and regardless of its microscopic degrees of freedom and even if, like the theory in Ref. [23], it 
escapes the instability argument of [12] . However, even if such an instability does exist, it may 
be avoided if w < — 1 for a sufficiently limited duration [21]. 
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